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Abstract. We give a non-perturbative construction of the fermionic projector in 
Minkowski space coupled to a time-dependent external potential which is smooth 
and decays faster than quadratically for large times. The weak and strong mass 
oscillation properties are proven. We show that the integral kernel of the fermionic 
projector is of Hadamard form, provided that the time integral of the spatial sup- 
norm of the potential satisfies a suitable bound. This gives rise to an algebraic 
quantum field theory of Dirac fields in an external potential with a distinguished 
pure quasi-free Hadamard state. 
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1. Introduction 

Hadamard states play an important role in quantum field theory because they are 
needed as a starting point for the perturbative treatment. Namely, interactions can be 
dealt with within the framework of perturbation theory only after a suitable regular¬ 
ization, yielding the so-called Wick polynomials. By using Hadamard states, one can 
devise a specific regularization scheme which is covariant and, thus, especially suited 
to be applicable also on curved backgrounds (see for example [mil! or the recent text 
book l^i. 

While the existence of Hadamard states can be shown abstractly with glueing con¬ 
structions (see [2Z1126I), these methods are not explicit. Therefore, the construction 
of Hadamard states remains an important task. In recent years, two methods for 
constructing Hadamard states have been developed: the pseudo-differential operator 
approach by Gerard and Wrochna [29l [301 EE] and the holographic method by Dappi- 
aggi, Moretti and Pinamonti laiSlITlEI. The first method gives a whole class of states. 
The second method does give a distinguished state, but it has the shortcoming that it 
relies on conformal methods which only apply to massless fields. 

Here we propose a new method for constructing Hadamard states for massive Dirac 
helds, which gives rise to a distinguished Hadamard state. The method is based on the 
recent functional analytic construction in [221123] of the fermionic projector in globally 
hyperbolic space-times. In the present paper, we show that the integral kernel of 
the fermionic projector yields a pure quasi-free Hadamard state. In order to make 
an explicit analysis feasible, we restrict attention to the setting of a time-dependent 
external potential in Minkowski space. We consider this work as a starting point for 
more general constructions in curved space-times, to which our methods should also 
apply. We point out that our method does not require any symmetries of space-times. 
However, if there are space-time symmetries, our state has these symmetries as well. 

More specifically, we show that the construction in inhnite lifetime in [23] applies 
to the Dirac equation in Minkowski space in the presence of an external potential, 
provided that the potential is smooth and decays suitably for large times. The main 
technical step is to prove that the so-called mass oscillation property holds. Assuming 
in addition a bound on the time integral of the spatial sup-norm of the potential, we 
show that the resulting fermionic projector is of Hadamard form (for an introduction 
to the Hadamard form see [33l El])- These results put the previous perturbative 
treatment of the fermionic projector in [nmanaiiHiisi] (see also the textbook [T5] i 
on a rigorous functional analytic basis. In particular, our results show that the nonlocal 
low and high energy contributions as introduced in m by a formal power series are 
indeed well-defined and smooth. 

In the remainder of the introduction, we state our results and put them into the 
context of the fermionic projector and of algebraic quantum held theory. 

The Mass Oscillation Property. In Minkowski space without an external poten¬ 
tial, the Dirac equation has plane-wave solutions. The sign of the frequency of these 
plane-wave solutions gives a splitting of the solution space into two subspaces, usually 
referred to as the positive and negative energy subspaces. This frequency splitting is 
important for the physical interpretation of the Dirac equation and for the construc¬ 
tion of a corresponding quantum held theory. Namely, choosing the vacuum state in 
agreement with the frequency splitting (Dirac sea vacuum), it is possible to reinterpret 
the negative-energy solutions in terms of anti-particle states. The plane-wave solutions 
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of positive and negative frequencies are then identified with creation and annihilation 
operators, respectively, which by acting on the vacuum state generate the whole Fock 
space. 

The above frequency splitting can still be used in static space-times (i.e. if a time- 
like Killing field is present). However, in generic space-times or in the presence of a 
time-dependent external potential, one does not have a natural frequency splitting. A 
common interpretation of this fact is that there is no distinguished ground state, and 
that the notion particles and anti-particles depend on the observer. Nonetheless, the 
construction of the fermionic projector as carried out non-perturbatively in |22l I23| 
does give rise to a canonical splitting of the solution space of the Dirac equation into 
two subspaces even in generic space-times. This also suggests that, mimicking the con¬ 
struction for the usual frequency splitting, there should be a canonical ground state of 
the corresponding quantum field theory, even without assuming a Killing symmetry. 
One of the goals of this paper is to construct this distinguished ground state in the 
presence of a time-dependent external potential in Minkowski space. 

We now recall a few basic constructions and definitions from |23) , always restricting 
attention to subsets of Minkowski space and to the Dirac equation 

{i(^ + 'B{t,x) - m)ijjrn{t,x) = 0 (1.1) 


in the presence of a smooth external potential S. Here m is the rest mass, and for 
clarity we add it as an index to the wave function. The construction differs considerably 
in the cases when the space-time has finite or infinite lifetime. A typical example of 
a space-time of finite lifetime is an open subset D of Minkowski space contained in a 
strip (—T, T) X for some constant T > 0 such that the surface {0} x is a Cauchy 
surface (for a general treatment of space-times of finite lifetime see [22]). In this case. 


one considers on the solution space of the Dirac equation (1.1) the usual scalar product 
obtained by integrating over the Cauchy surfac^ 


— 27r / fim'^\[t=0,x) d X 


( 1 . 2 ) 


as well as the space-time inner product 

x d X (f’3) 

Jn 

(here -<fi\4>>- is the spin scalar product, which is often denoted by with the adjoint 
spinor where the dagger means complex conjugation and transposition). 

Finite lifetime implies that the space-time inner product is bounded in the sense that 
there is a constant c > 0 such that 


\<4>m\'fim>\ < C 


for all Dirac solutions (and 


\\m llV’m||m (1-^) 

is the norm corresponding to the scalar 
product (1.2)). This in turn makes it possible to represent the space-time inner product 
in terms of the fermionic signature operators, meaning that there is a unique bounded 
symmetric operator § such that 

— {fim \ ^'fira)m 


^The factor 27r might seem unconventional. This convention was first adopted in m to simplify 
some formulas. 
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(here the tilde indicates that an external potential ® is present, whereas the corre¬ 
sponding objects in the Minkowski vacuum are denoted without a tilde). Then the 
positive and negative spectral subspaces of S give rise to the desired splitting of the 
solution space. 

The above construction fails in space-times of infinite lifetime because the time 
integral in (1.3) will in general diverge. The way out is to consider families of so¬ 
lutions of the family of Dirac equations (1.1) with the mass parameter m 


varying in an open interval I. We need to assume that I does not contain the origin, 
because our methods for dealing with inhnite lifetime do not apply in the massless 
case m = 0 (this seems no physical restriction because all known fermions in nature 
have a non-zero rest mass). By symmetry, it suffices to consider positive masses. Thus 
we choose 

I := (mi, mu) C M with parameters mL,mR > 0 . (1-5) 

We always choose the family of solutions (V’m)me/ ™ the class x I,SJi) 

of smooth solutions with spatially compact support in Minkowski space M which 
depend smoothly on m and vanish identically for m outside a compact subset of I. 


On such families of solutions, we can impose conditions analogous to (1.4) by suitably 


integrating over m. We here give the condition which is most relevant for applications 
(for a weaker version, which will also arise in intermediate steps of our proofs, see 


Definition 2.3 below). 


Definition 1.1. The Dirac operator il^ +'B has the strong mass oscillation prop¬ 
erty in the interval I (see (|1.5 [) ) if there is a constant c > 0 such that 


< / firn dm I / ifrn' dm'> 


< C 


dm 


( 1 . 6 ) 


for all families of solutions {ipm,)mei, {4>m)mei £ x fiSM). 


The point is that we integrate over the mass parameter before taking the space-time 
inner product. Intuitively speaking, integrating over the mass parameter generates a 
decay of the wave function, making sure that the time integral converges. 

As shown in |23l Section 4], the strong mass oscillation property gives rise to the 
representation 

<j4>mdm\ J firn' dm'> = J {4>jn \ §m V’m)m dm , 
which for every m ^ I uniquely defines the fermionic signature operator S^n- This 


operator is bounded and symmetric with respect to the scalar product (1.2). Moreover, 
it does not depend on the choice of the interval I. Now the positive and negative 
spectral subspaces of the operator again yield the desired splitting of the solution 
space. 


The remaining crucial question is whether the inequality (1.6) holds in the presence 

Theorem 5.1] 


of an external potential. In the case of a static potential, it is shown in 
that the mass oscillation property holds, and that the fermionic projector gives back 
the usual frequency splitting. In particular, if 23 vanishes, the fermionic projector 
gives the usual vacuum ground state. With these results in mind, we here consider the 
time-dependent situation. Also in this setting, we give an affirmative answer to the 
above question, provided that the potential has suitable decay properties at infinity: 
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Theorem 1.2. Assume that the external potential 3 is smooth and for large times 
decays faster than quadratically in the sense that 

mt)\c^ < (1.7) 

for suitable constants e,c > 0. Then the strong mass oscillation property holds. 


The C^-norm in (1.7) is defined as follows. We denote spatial derivatives by V and use 
the notation with multi-indices, i.e. for a multi-index a = (ai,..., ap) we set V“ = 




and denote the length of the multi-index by |q:| = p. Then the spatial C^^-norms 


of the potential are defined by 


\3{t)\Qk := max sup 
l"l<^ xeiR3 




( 1 . 8 ) 


where | . | is the sup-norm corresponding to the norm |(/>p := (!)>- on the spinors. 


The Fermionic Projector and the Hadamard Form. Using the result of the 
previous theorem, the fermionic projector P is defined for a fixed mass parameter m ^ I 
by (see [22l Dehnition 3.7] and [23l Definition 4.5]) 

P ~X{—oo,0){^m) km : (1-9) 

where X(_oo,o) (Sm) is the projection onto the negative spectral subspace of the fermionic 
signature operator, and km is the causal fundamental solution (for basic dehnitions 
see Section [2.1| below). The fermionic projector can be represented as a bi-distribu- 
tion P{x,y) with x,y G M (see [22l Section 3.5] and [23l Section 4.3]), which satishes 
the Dirac equation and is symmetric, i.e. 

{i0^ + 3{x) — m) P{x,y) = 0 (1-10) 

P{x,yr = P{y,x) (1.11) 


(where P{x,y)* = 7 ^P(x,y)l 7 ^ is the adjoint with respect to the spin scalar prod¬ 
uct 7 . 1 . 7 ). 

Knowing the singularity structure of the bi-distribution P{x,y) is important for ap¬ 
plications (point-splitting method, Wick polynomials, renormalization, etc.). There¬ 
fore, we shall establish that the fermionic projector is of Hadamard form. In our 
setting, this is tantamount to proving that the bi-distribution P{x, y) is of the form 
(see [31] or [32l page 156]) 

P(x,y) = lim i^ (^^-^Y^ + V{x,y) log a e{x,y) + W{x,y'^ , (1.12) 

e\o ^ \ae[x,y) j 


where 

y) := (y - xy {y - x)j - ie {y - x)° , (1.13) 

and U, V and W are smooth functions on M x M taking values in the 4 x 4-matrices 
acting on the spinors (we always denote space-time indices by latin letters running 
from 0,..., 3). For clarity, we point out that on a manifold, the function can be 
defined locally in a geodesically convex neighborhood of a point x G M, making it 
necessary to distinguish between the local Hadamard form (i.e. a local representation 
of the form (1.12)) and the global Hadamard form (implying that the singularities 
in (1.12) are the only singularities of the bi-distribution). For these subtle issues, we 
refer the reader to [all [271 [36]. In our setting of Minkowski space, there is one global 
chart, and the distance function is defined globally by (1.13). For this reason, we 
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do not need to make a distinction between the local Hadamard form and the global 
Hadamard form. 

In space-times of finite lifetime, in general the fermionic projector is not of Hadamard 
form. The first counter examples were constructed in |10j . where it is shown that in 
ultrastatic space-times of finite lifetime, the Hadamard condition is in general violated. 
Other counter examples are so-called simple domains as introduced in m Sections 2.1 
and 1.4]. In such simple domains, the fermionic projector in the massless case is an 
operator of finite rank with a pointwise bounded integral kernel, clearly not being of 
Hadamard form. 

In space-times of infinite lifetime, the situation is better at least in the ultrastatic 
case. Namely, in |23[ Section 5] it is shown that the fermionic projector in ultra¬ 
static space-times is composed of all negative-energy solutions of the Dirac equation. 
Therefore, P{x, y) coincides with the bi-distribution constructed from the frequency 
splitting, which in |32] was shown to be of Hadamard form. Apart from this specihc 
result, it is unknown whether or in which space-times the fermionic projector is of 
Hadamard form. 

The main result of this paper is to show that in a time-dependent external potential 
in Minkowski space, the fermionic projector is indeed of Hadamard form, provided 
that the potential is not too large: 


Theorem 1.3. Assume that the external potential “B is smooth, and that its time 
derivatives deeay at infinity in the sense that (1.7) holds and in addition that 


f 


\dlB{t)\(jo dt < oo for all p G N 


(with the C^-norm as defined in (1.8)J. Moreover, assume that the potential satisfies 
the bound 


/ \‘B{t)\co dt < V2 - I . 

J —OO 

Then the fermionic projector P{x, y) is of Hadamard form. 


(1.14) 


We note that the property of the bi-distribution P{x,y) to be of Hadamard form 
can also be expressed in terms of the wave front set (see for example [Ml ED]). Also, 
the smooth functions in drill ) can be expanded in powers of the Minkowski dis¬ 
tance a{x,y) = {y- xfi {y -x)j, 

OO OO OO 

U{x,y) = ^Un{x,y)a'^ , V{x,y) = '^Vn{x,y) , W{x,y) = ^Wn{x,y) . 

n=0 n=0 n=0 


The coefficients of this so-called Hadamard expansion can be computed iteratively using 
the method of integration along characteristics (see [MIES] or 121). In Minkowski space, 
the light-cone expansion [MUM] gives a systematic procedure for computing an infinite 
number of Hadamard coefficients in one step. This procedure also makes it possible 
to compute the smooth contributions to P{x,y), giving a connection to fermionic loop 
corrections in quantum held theory (see m § 8.2 and Appendix D]). 

We finally remark that, introducing an ultraviolet regularization, the fermionic pro¬ 
jector gives rise to a corresponding causal fermion system (for details see [22l Sec¬ 
tion 4]). In this context, the result of Theorem 1.3 gives a justification for the for¬ 
malism of the continuum limit as used in mm for the analysis of the causal action 
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principle. We also refer the interested reader to the introduction to causal fermion 
systems [20] , 


Quantum Fields and Hadamard States. Using the standard notation in quantum 
field theory, the objective of the quantization of the Dirac field is to construct field 
operators and '!'(?/)* acting on a Fock space dfpock together with a suitable ground 
state |0). The field operators should satisfy the canonical anti-commutation relations 


{T“(x),T^(y)} = {T“(x)*,T^(yr} = 0 


(1.15) 

(1.16) 


where the Greek indices running from 1,..., 4 denote Dirac spinor indices (we always 
work in natural units h = c = 1). For clarity, we point out that, due to the conjugation 
of the operator ^^{y) in ( 1.15[ ), the spinor index /3 should be thought of as the index of 
a dual spinor; this is why it appears on the right side as lower index. In the presence of 
a time-dependent external potential, there is no distinguished ground state. But there 
is common agreement that for a physically sensible theory the ground state |0) should 
be chosen such that the two-point function (0| 'I'(x) T(y)* |0) is a bi-distribution of Ha¬ 
damard form. In this paper, we shall achieve this goal by arranging that the two-point 
function coincides with the bi-distribution of the fermionic projector, i.e. 

(0|T“(x)T^(y)*|0) = -(P(x,y));. 

In order to give the above formulas a mathematical meaning, one needs to “smear 
out” the field operators and work with operator-valued distributions (see for exam¬ 
ple [39l El El 0] ). Formally, this is accomplished by setting 

T( 5 ) = [ T(x)“ 5 r(x)„ and ^*(/) = [ (T(x)")*/(a:)" , 

Jjt Jji 

where g and / are smooth and compactly supported co-spinors and spinors, respec¬ 
tively. We do not aim at defining the pointwise field operators, but instead we work 
exclusively with the smeared field operators 4'(g() and 4'*(/) (for basic definitions see 
Section]^. Moreover, instead of considering vacuum expectation values, in the alge¬ 
braic formulation of quantum field theory one prefers to work with a pure quasi-free 
state cn, making it unnecessary choose a representation of the field algebra on the Fock 
space. Given a state w, a corresponding representation of the field algebra is obtained 
by applying the GNS construction, also making it possible to recover w as a vacuum 
expectation value. A quasi-free state for which the two-point function is of Hadamard 


form (1.12) is called a Hadamard state. 


Using the algebraic language, we prove the following result: 


Theorem 1.4. There is an algebra of smeared fields generated by ^{g), 4'*(/) together 
with a pure quasi-free state u with the following properties: 

(a) The canonical anti-commutation relations hold: 

{^{g), T*(/)} = <g* \kmf>, {'k(ff), 'k(5')} = 0 = {T*(/), T*(/')} . (1.17) 

(b) The two-point function of the state is given by 

w(T(ff)T*(/)) =- If g{x)P{x,y)f{y)d‘^xd'^y . 

JJ MxM 
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The main step in the proof is to use the spectral projection operators X(-oo,0)('5m) 
and X[o,oo)('5m) to construct a positive operator R, making it possible to apply Araki’s 
results in [T] to obtain the desired quasi-free state. 

We finally put our result into the context of other methods for constructing Hada- 
mard states. First, there is the method of glueing the physical space-time to an 
ultrastatic space-time and using that the Hadamard property is preserved under time 
evolution (see 1271 [26]). This method shows the existence of Hadamard states in 
every globally hyperbolic space-time and gives a constructive procedure for a class of 
Hadamard states. Another method is to work with pseudo-differential operators [29], 
again giving a whole class of Hadamard states. A method which distinguishes one 
specific Hadamard state using asymptotic symmetries at null infinity is given in |^. 
Our method gives a unique distinguished Hadamard state even in the generic time- 
dependent setting in Minkowski space. Moreover, this method is constructive in the 
sense that the bi-distribution P{x, y) and its Hadamard expansion can be computed 
explicitly (see IMHHl). Our results exemplify that the construction of the fermionic 
projector in [23] is a promising method for constructing a distinguished Hadamard 
state without any symmetry assumptions, hopefully even in generic globally hyperbolic 
space-times. 


2. Preliminaries 


2.1. Dirac Green’s Functions and the Time Evolution Operator. Let Ji be 

Minkowski space, a four-dimensional real vector space endowed with an inner product 

of signature (-1-). The Dirac equation in the Minkowski vacuum (i.e. without 

external potential) reads 

{i^ — m) Jp{x) = 0 , 

where we use the slash notation with the Feynman dagger 0 := We always work 
with the Dirac matrices in the Dirac representation 


0 = r ® 

V 0 -lc2 



a 

0 


(and a are the three Pauli matrices). The wave functions at a space-time point x take 
values in the spinor space Sx, a four-dimensional complex vector space endowed with 
an indefinite scalar product of signature (2,2), which we call spin scalar product and 
denote by 


4 

^ Sa V’"(a^)V"(a;) , Si = S2 = 1, S3 = S4 = -1 , 

Q=1 


where is the complex conjugate wave function (this scalar product is often written 
as 'ipcf) with the so-called adjoint spinor We denote the space of smooth wave 

functions by C°° {M , SJ/t ), whereas {M , SJi) denotes the smooth and compactly 
supported wave functions (here SJi is the spinor bundle over Minkowski space with 
hbers SxJi)- On the spaces of wave functions, one can introduce a Lorentz-invariant 
pairing by integrating the spin scalar product over space-time. 


<.|.> : C°^{Jl,SJl) xC^{Ji,SJl) ^C, 







( 2 . 1 ) 
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In what follows, the mass parameter of the Dirac equation m will not be fixed. It 
can vary in an open interval I := {mL^mn) with mi,mu > 0. In order to make this 
dependence explicit, we often add the mass as an index. Moreover, we consider an 
external potential “B, which we assume to be symmetric with respect to the spin scalar 
product. 


^|J I 4>>~x = I B(x) for ail X £ Ji and G Sx 

Then the Dirac equation becomes 


( 2 . 2 ) 


(2? — m) 'ijjm = 0 with V := ilj) + B . 


(2.3) 

Since the Dirac equation is linear and hyperbolic (meaning that it can be rewrit¬ 
ten as a symmetric hyperbolic system), its Cauchy problem for smooth initial data 
is well-posed, giving rise to global smooth solutions. Moreover, due to finite prop¬ 
agation speed, starting with compactly supported initial data, we obtain solutions 
which are spatially compact at any time. We denote the space of such smooth wave 
functions with spatially compact support by S.M). Using the symmetry as¬ 
sumption (2.2), for any solutions G SM) of the Dirac equation the 


vector field (t>my is divergence-free; this is referred to as current conservation. 

Applying Gauss’ divergence theorem, this implies that the spatial integral 


(V’mII^ •— 27r / 4^m'^\(t,x) d X 


(2.4) 


is independent of the choice of the space-like hypersurface labelled by the time pa¬ 
rameter t. This integral defines a scalar product on the solution space corresponding 
to the mass m. Forming the completion, we obtain a Hilbert space, which we denote 
by {Bim, (-I-)™)- The norm on Bm is ||.||m- 

The retarded and advanced Green’s operators s^ and s)^ are mappings (for details 
see for example m) 


Their difference is the so-called causal fundamental solution km, 

km ^ (s]^ - s^) : C^{Ji,SJl)^C^iM,SM)nBm. (2.5) 

ZTTl 

These operators can be represented as integral operators with a distributional kernel, 
for example, 

{km(p)ix)= kmix,y)(j){y) d'^y . 

Jji 

Leaving out the tilde always refers to the special case 23 = 0. 

The operator km can be used for constructing a solution of the Cauchy problem. To 
this end, we always work in the foliation A) = {(t, x) | x G of constant time Cauchy 
hypersurfaces in a fixed reference frame (t, x). For clarity, we denote the Hilbert space 
of square integrable spinors at time t with the scalar product (2.4) by (FCt, (.|.)|t). 
Moreover, we denote a wave function ■0 at time thy 'ip\t (we use this notation both for 
the restriction of a wave function in space-time and for a function defined only on the 
hyperplane jV*). 


Proposition 2.1. The solution of the Cauchy problem 
(P - m) V’m = 0 , = lAo e 




( 2 . 6 ) 
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has the representation 

V’r 


For the proof see for example |22l Section 2]. 

Moreover, the operator km can be regarded as the signature operator of the inner 


i{x) = 2tt km{x,ito,y))-f°-ipo{y)(fy. 


product (2.1) when expressed in terms of the scalar product (2.4). 


Proposition 2.2. For any il^m £ O'nd (j) G SJi), 

{'4’m I km(f)m — ■ (2-7) 

For the proof we refer to [9l Proposition 2.2] or [22l Section 3.1]. 

The unique solvability of the Cauchy problem allows us to introduce the group of 
time evolution operators as follows. According to Proposition 2.1, for given initial 


data ipQ G the Cauchy problem (2.6) has a unique solution if^m £ 

SM) n Jim- Evaluating this solution at some other time t gives a map- 


ping Um° '■ '00 0m|i- Since the scalar product (2.4) is time independent, the 
time evolution operator Um° is isometric. Thus by continuity, it extends uniquely to 
an isometry 


U, 


t,to 


Ji, 


to 


Jit 


Since to can be chosen arbitrarily and the Cauchy problem can be solved forward and 
backward in time, this isometry is even a unitary operator. Moreover, these operators 
are a representation of the group (M, meaning that 


= 1 


and 


jjt^kt' Tjt\t _ 

m m 


Proposition 

kernel, 


2.1 




immediately yields the following representation of Um with integral 


010 (0) = / ^m*(0 x) 0|t(f) (fx , 

4r3 

Um\y,x) = 2-71 km{{t',y), it,x)) 7 ° . 


( 2 . 8 ) 

(2.9) 


2.2. The Mass Oscillation Property. We denote the families of smooth wave 
functions with spatially compact support, which are also compactly supported in I, 
by X I,SM). The following construction shows that within this class, there 

are families {'ipm)m£i such that for every m ^ I, the wave function 0^ is a solution of 
the Dirac equation (1.1): Let 0° G C^Q{J{to x I, SM) be a family of smooth and com¬ 
pactly supported functions on jV x / (for example of the form 0^(x,m) = r]{m)x{x) 
with T] G C'q“(/) and x G C“(7V, SJi)). Solving for every m £ I the Cauchy problem 

{i^+'B- m)'ll)m = 0 , 0mU = 0m 

(where again 0^(x) = 0‘^(x, m)), we obtain a family {'ipm)m,ei of solutions of the Dirac 
equation for a variable mass parameter m G / in the desired class C^q^JI x I, SJi). 

The space of families of Dirac solutions in the class C^q{JI x I, SJi) are denoted 
by TC°°. On Ji°° we introduce the scalar product 


(010) = / (0m|0m)m dm , 


(2.10) 
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where dm is the Lebesgue measure (and V’ = (V'm)me/ and cf) = {4’m)mei are families 
of Dirac solutions for a variable mass parameter). Forming the completion yields the 
Hilbert space (IK, (.|.)) with norm ||.||. Then K°° can be regarded as the subspace 

= C^q{JI X I,SJl)nJi. ( 2 . 11 ) 

On K, we introduce the operator of multiplication by m, 

T : K-^K, {Tip)m = mi/jm ■ 

Obviously, this operator preserves the support properties, and thus 

Til , n_roo , rirOO 

1 ^ ^ 

Moreover, it is a symmetric operator, and it is bounded because the interval I is, i.e. 

T* = T £ L(K) . 

Integration of ^pm. over m gives another operator 

p : ^ SJl) , p^jJ = J il^mdm. (2-12) 

We point out for clarity that pi/i no longer satisfies a Dirac equation. The following 
notions were introduced in [23], and we refer the reader to this paper for more details. 

Definition 2.3. The Dirac operator D = i^ + T> on Minkowski space Ji has the weak 
mass oscillation property in the interval I = {mL,mR) with domain K°° if the 
following conditions hold: 

(a) For every £ K°°, the function is integrable on M. Moreover, 

there is a constant c = c{'ifi) such that 

|<p'i/i|p(?!)>| < c ||(/>|| for all (/) G K°° . (2.13) 

(b) For all if,4>£ K“, 

<prV’|p<?i> = <pi/>|pT(/)> . (2-14) 

Definition 2.4. The Dirac operator D = il^+'B on Minkowski space Ji has the strong 
mass oscillation property in the interval I = {mL,mji) with domain K°° if there 
is a constant c > 0 such that 

|<p'0|p(/)>| < 0^ \\4)m,\\m\\'ipm\\mdm for all i/i, (/) G K“ . (2.15) 

The following theorem is proved in [231 Theorem 4.2, Proposition 4.3 and Theorem 4.7]. 

Theorem 2.5. Assume that the Dirac operator D has the strong mass oscillation 
property in the interval I = (m^,, m/j). Then there exists a family of linear operators 
(Sm)meJ with §m £ L(Km) which are uniformly hounded, 

sup ||Sm|| < OO , 

m&I 

such that 

<pf;\p(j)> = J{i:m\§m4>m)mdm for all £'K°° . (2-16) 

The operator §rn is uniquely determined for every m £ I by demanding that for 
all £ IK°°, the functions {tpm\^m4>m)m o-re continuous in m. Moreover, the 
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operator Sm is the same for all choices of I containing m. Finally, there is a bi¬ 
distribution V G X J/i) such that the operator P defined by 

p ■■= -X(-ocfl){§m)~km : (2.17) 

has the representation 

= (2.18) 
(where cf) = is the usual adjoint spinor). 


The operator P is referred to as the fermionic projector. We also use the standard 
notation with an integral kernel P{x,y), 


<(j)\Pi»= P4>{x)\P{x,y)i>{y)yxd xd y 

JJ jHxM 

{P'if){x)= ! p{x,y)'if{y)d^y , 

Jji 


where P{. ,.) coincides with the distribution V in (2.18). 

2.3. The Lippmann-Schwinger Equation. The Dirac dynamics can be rewritten 
in terms of a symmetric operator H. To this end, we multiply the Dirac equation (2.3) 
by 7 ® and bring the t-derivative separately on one side of the equation, 

idtijm = H'f’m , where H := — 7 °(i 7 V + T> — m) (2-19) 

(note that 'y^dj = 7 ^ 9 * + 7 V). We refer to (2.19) as the Dirac equation in Hamiltonian 
form. The fact that the scalar product (2.4) is time independent implies that for any 
two solutions SJi) n !Km, 

0 = dt{(t>rn I V^m)m — i{{H 4 *m \ {'Pm \ H''Pm)m) j 

showing that the Hamiltonian is a symmetric operator on 34^- The Lippmann- 
Schwinger equation can be used to compare the dynamics in the Minkowski vacuum 
with the dynamics in the presence of an external potential. We denote the time evo¬ 
lution operator in the Minkowski vacuum by Um°. 

Proposition 2.6. The Cauchy problem (2.6) has a solution ifm which satisfies the 
equation 


P’mlt = + i u(yf (7°T lAm) 1^ dr , 

Jto 

referred to as the Lippmann-Schwinger equation. 


( 2 . 20 ) 


Proof. Obviously, the wave function 'tpmlt given by (2.20) has the correct initial condi¬ 
tion at t = to- Thus it remains to show that iprnlt satisfies the Dirac equation. To this 
end, we rewrite the Dirac equation in the Hamiltonian form (2.19), and separate the 
vacuum Hamiltonian H from the term involving the external potential, 

{idt — H) ifm = —7^'^ 'Pm with H = —i^^pV +. ( 2 - 21 ) 

Applying the operator idt—H to (2.20), and observing that the time evolution operator 
maps to solutions of the vacuum Dirac equation, only the derivative of the upper limit 
of integration contributes, 

{idt - H) prnV = ( 7 °® pra) = - 7 °'® Pm\t , 


so that ( 2 . 21 ) is indeed satisfied. 


□ 
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3. The Mass Oscillation Property in the Minkowski Vacuum 


Since Minkowski space is ultrastatic, it is known from |23l Section 5] that the Dirac 
operator i(^ satisfies the weak and strong mass oscillation properties. Moreover, the 
decomposition of the solution space into the positive and negative spectral subspaces 
of the fermionic signature operator reduces to the usual frequency splitting (see |23l 
Theorem 5.1]). We now reproduce these results giving more explicit proofs. These 
explicit results and formulas will be essential for the subsequent treatment of time- 
dependent external potentials in Section 

Basically, the mass oscillation property in the Minkowski vacuum can be proved 
easily using Fourier methods. Here we shall give two different approaches in detail. 
The method of the first proof (Section |3.1[ ) is instructive because it gives an intuitive 
understanding of “mass oscillations”. However, this method only yields the weak mass 
oscillation property. The second proof (Section |3.2[ ) is more abstract but also gives the 
strong mass oscillation property. 


3.1. Proof of the Weak Mass Oscillation Property using Mass Derivatives. 

We again consider the foliation Jit = {(^j jl) | ^ £ 1^^} of constant time Cauchy hy¬ 
persurfaces in a fixed reference frame {t, x) and a variable mass parameter m in the 
interval I = with mi, mu > 0. The families of solutions if: = of 

the Dirac equations {i(p — = 0 are contained in the Hilbert space (IK, (.|.)) with 

the scalar product (2.10). Moreover, the subspace 1K°° C IK is given by ( 2.11| ). 

For what follows, it is convenient to work with the Fourier transform in space, i.e. 

f d^x , 'il;{t,x)= [ 

Jr 3 JR3 (27r)'^ 

Then a family of solutions ijj G 1K°^ has the representation 

^m{t, k) = c+{k,m) + c-{k, m) 


i)it,k) = 


k) e' 


ikx 


for all m G / 


(3.1) 

with suitable spinor-valued coefficients c±(fc, m) and u}{k, m) := \J\k\‘^ + m?- Integrat¬ 
ing over the mass parameter, we obtain a superposition of waves oscillating at different 
frequencies. Intuitively speaking, this leads to destructive interference for large t, giv¬ 
ing rise to decay in time. This picture can be made precise using integration by parts 
in m, as we now explain. Integrating (3.1) over the mass and applying the operator p, 

"*‘^*) dm 


(2.12), we obtain 

P'4’it,k) = 


(c+ + c_ e“ 


11 t dm^ 
i 


t 


d. 


(c+ dm(^ 

C+ 
dmU^ 


- C_ 


) dm 


—iujt 


- d„ 


C- 

dm,X} 


Aojt 


dm 


(we do not get boundary terms because G 'K°° has compact support in m) 
dmUJ = m/ui, we conclude that 


With 


P^{t,k) = - 


t 


’I 


dr, 


m / 


U) C- 


m 


Aujt 


dm . 


Since the coefficients c± depend smoothly on m, the resulting integrand is bounded 


uniformly in time, giving a decay at least like 1/t, i.e. 


< 


1/t. Iterating 
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this procedure, one even can prove decay rates < 1/t^, 1/t^,... The price one pays is 
that higher and higher powers in uj come up in the integrand, which means that in 
order for the spatial Fourier integral to exist, one needs a faster decay of c± in \k\. 
Expressed in terms of the initial data, this means that every factor \/t gives rise to an 
additional spatial derivative acting on the initial data. This motivates the following 
basic estimate. 


Lemma 3.1. For any iIj G !K°°, there is a constant C = C{mL) such that 

2 


|(pV')k|L < 


Cl 


1 +t2 




m£l 


'mjU=0||^2,2 ? 


(3.2) 


b=0 


where ||.||t is the norm corresponding to the scalar product 

(.|.)|t:=27r f SM) x SJi) 


(which is similar to (2.4), but now applied to wave functions which do not need to be 
solutions), and ||.11^4/2,2 is the spatial Sobolev norm 


|2 

I VF2'2 




(3.3) 


a with |a| < 2 ^ 

where a is a multi-index. 

The absolute value in (3.3) is the norm | . | := on the spinors. If we again 

identify all spinor spaces in the Dirac representation with C^, this simply is the stan¬ 
dard Euclidean norm on C^. 

The proof of this lemma will be given later in this section. Before, we infer the weak 
mass oscillation property. 

Corollary 3.2. The vacuum Dirac operator i(^ in Minkowski space has the weak mass 
oscillation property with domain (2.11). 

Proof. Eor every fjjf) € the Schwarz inequality gives 


I<pV'Ip</»I = ^ 


((pV’)k |7° {p(()U)tdt 


< 


AlipMlf dt. (3.4) 


Applying Lemma |3.1| together with the estimate 


(P<?5’)|i||i = // {4>m\t\4>m'\t)^dmdm' 


Ixl 

1 

< - 

2 .1.1 ixi 


'm\\m 


+ \\4>m'\\m') dm dm' = \I\ 


we obtain inequality (2.13) with 

2 


c = C\I \2 sup ^ ||5^(V'm)|t=o||ry2,2 


m£l 


b=0 


'_oo 1 + 


dt < oo . 


(3.5) 


The identity (2.14) follows by integrating the Dirac operator by parts, 

<pTil;\p(j)> = <pVij;\p(j)> = <Vp'ij;\p(j)> = / ^Vp'i^pfiyx d^x 

Jji 


(*) 


(3.6) 


'Ji 


~<p'il)\Vp(l))^x d^x = <pilj\Vp(l)> = <p-i/)|pT(/)> . 
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In (*), we used that the Dirac operator is formally self-adjoint with respect to the 
inner product <.|.>. Moreover, we do not get boundary terms because of the time 
decay in Lemma 3.1 □ 


The remainder of this section is devoted to the proof of Lemma 3.1 Specializing the 
result of Proposition 2.1 to the Minkowski vacuum, we can express the solution i/’m 
of the Cauchy problem in terms of the causal fundamental solution km- In order to 
bring km into a more explicit form, we use (2.5) together with formulas for the advanced 
and retarded Green’s functions. Indeed, these Green’s functions are the multiplication 
operators in momentum space 


Sm{k) = lini 


-|- m 


e\o k"^ — m? — iek^ 


and Sm{k) = lim 


: -|- m 


s\o k"^ — m? + iek^ 


(with the limit e 0 taken in the distributional sense, and where the vector k is the 
four-momentum). We thus obtain in momentum space 


km{p) = 


27ri 


I -|- m) lim 
£\0 


1 


=- ii) + m) lim 

27ri ^ e\o 


p2 _ ^2 _ jgpO 


p2 _ jy^2 _ jg 


p2 _ ^2 _|_ j^^pO 

1 


p2 _ jy^2 _|_ 


e{p^) 


(where e denotes the step function, and for notational clarity we denoted the momen¬ 
tum variables by p). Employing the distributional equation 


lim ( — -^ 

e\o \x — ie X + ie 

we obtain the simple formula 


= 27ri 5{x) , 


km{p) = {p + m) 5{p^ - n?) e(p°) 


(3.7) 


It is convenient to transform spatial coordinates of the time evolution operator to 
momentum space. First, in the Minkowski vacuum, the time evolution operator can 
be represented as in (2.8) with an integral kernel (y, x) which depends only on the 
difference vector y — x. We set 

U*’^'{k):= f d^y . 


Combining (2.9) with (3.7) yields 

roo 

+ m) 7 ° 5{k^ - e(w) duj 


U^^^'ik) = / 

J —C 

Carrying out the w-integral, we get 


where we set 


[/‘’*'(fe) = ^n±(fe) 
± 


n±(fe) :=± 


2uj{k) 




+ m) 7° 


(3.8) 


(3.9) 


with ijj{k) = y \k\^ + m? and = (±a;(/c),/c) . 
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Moreover, applying Plancherel’s theorem, the scalar product (2.4) can be written in 
momentum space as 




m \ Ymjm 


(2vr)-2 f 

Jrs 


k) I 7° k)>- (fk . 


The unitarity of the time evolution operator in position space implies that the ma¬ 
trix (k) is unitary (with respect to the scalar product (., .)c 2 = -<. . )^), mean¬ 

ing that its eigenvalues are on the unit circle and the corresponding eigenspaces are 
orthogonal. It follows that the operators n±(A;) in (3.8) are the orthogonal projection 
operators to the eigenspaces corresponding to the eigenvalues \ i.e. 

70 n :70 = n, and n,(fc) ^ n.(fc) for s, s'G {+,-} . 


(these relations can also be verihed by straightforward computations using (3.9)). 

The next two lemmas involve derivatives with respect to the mass parameter m. 
For clarity, we again denote the m-dependence of the operators by the subscript m. 

Lemma 3.3. The time evolution operator in the vacuum satisfies the relation 


where 


(f - f') Uii'ik) = f^VilUk) + Wli'ik) , 

(3.10) 

Vpik) = E i (^± + 

(3.11) 

H7'W = E)(t^ j)■ 

(3.12) 

and Wm are estimated uniformly by 


l|i7"(fc)ii + i|iF7'WII<c(i + !h), 

V m) 

(3.13) 


where the constant C is independent of m, k, t and t' (and 
2 X 2-matrices). 


is any norm on the 


Proof. First, we generate the factor t — t' hy differentiating the exponential in (3.8) 
with respect to tv, 

(t - t') U‘J(k) = T.^±{k){±iP . 

Next, we want to rewrite the w-derivative as a derivative with respect to m. Taking 
the total differential of the dispersion relation — |A:p = for fixed fc, one finds 
that 

d io d 
du) m dm 

Hence 


(3.14) 


^ ^ ^ 


d 

m dm 




— n± '' — 


m 


E 


d 

±i — n± 

am L m 
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Computing the operators in the round brackets using (3.9) gives the identities (3.11) 
and (3.12). Estimating these formulas, one obtains bounds which are at most linear 
in |A;|, proving (3.13). □ 

This method can be iterated to generate more factors of t — In the next lemma, 
we prove at least quadratic decay in time. For later use, it is preferable to formulate 
the result in position space. 

Lemma 3.4. The time evolution operator in the vacuum has the representation 

= (t - t')2 

with operators 

{J{,, 

which are bounded uniformly in time by 

W^m {4>)\\t + ll-Bm {4>)\\t + \\C^ {4>)\\t < c ||</>||vi/2,2 , (3.16) 

where c is a constant which depends only on m. 


Proof. A straightforward computation using exactly the same methods as in Lemma 3.3 
yields the representation 






where the operators Am , Bm and Cm are bounded by 


|AV(t)ll + 


+ l|cV(t)||<-(i + S + lT 


m 




(3.17) 


(3.18) 


with a numerical constant C > 0. We remark that, compared to (3.10), the right 
of (3.18) involves an additional 1/m. This prefactor is necessary for dimensional 
reasons, because the additional factor t — t' in ( 3.17| ) (compared to (3.10)) brings in 
an additional dimension of length (and in natural units, the factor 1/m also has the 
dimension of length). The additional summand |/cp/m^ in (3.18) can be understood 
from the fact that applying (3.14) generates a factor of uj/m which for large \k\ scales 
like \k\/m. 

Translating this result to position space and keeping in mind that the vector k 
corresponds to the derivative —iV, we obtain the result. □ 


Proof of Lemma 3.1 


First of all, the Schwarz inequality gives 


(Pl/)|i||i < W'lfmWmdrn < y/\T\ ||l/|| . 


Thus it remains to show the decay for large t, i.e. 
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We apply Lemma 3.4 and integrate by parts in m to obtain 


1 


(pV')k = / V'mk=o dm = ^ I V’m|t=o dm 


II 

1 


J (^m (5mV'm|i=0) “ (9mV'm|i=o) + V’m|l=o) dm . 


Taking the norm and using (3.16) gives (3.19). 


□ 


We finally note that the previous estimates are not optimal for two reasons. First, 


the pointwise quadratic decay in (3.2) is more than what is needed for the convergence 


of the integral in (3.5). Second and more importantly, the Schwarz inequality (3.4) 


does not catch the optimal scaling behavior in k. This is the reason why the constant 
in (2.13) involves derivatives of i/jm (cf. ( |3.5[ )), making it impossible to prove inequal¬ 
ity (2.15) which arises in the strong mass oscillation property. In order to improve 


the estimates, one needs to use Fourier methods both in space and time, as will be 
explained in the next section. 

3.2. Proof of the Mass Oscillation Property using a Plancherel Method. 
Theorem 3.5. The vacuum Dirac operator in Minkowski space has the strong mass 


oscillation property with domain (2.11). 


Our proof relies on a Plancherel argument in space-time. It also provides an alternative 
method for establishing the weak mass oscillation property. 


Proof of Theorem 3^. Let if = be a family of solutions of the Dirac 

equation for a varying mass parameter in the Minkowski vacuum. Using Proposi¬ 


tion 2.1, one can express V'm in terms of its values at time t = 0 by 
ipm{x) = 2 tt kjnix, (0, y)) 7 ° V'm|t=o(y) d^y . 


We now take the Fourier transform, denoting the four-momentum by k. Using (3.7), 
we obtain 

fimik) = 2Trkm{k) 'j^'ip^ik) 

= 271 5{k^ — m^) e{k^) (Jf + m) 'y^'ijj^{k) , 

where denotes the spatial Fourier transform of V’m|t=o- Obviously, this is a 

distribution supported on the mass shell. In particular, it is not square integrable 


over 


Integrating over m, we obtain the following function 

{p'ip)ik) = 2 -k xi{m) ^ e{k°) + m) 

2m 




(3.20) 


where m now is a function of the momentum variables. Since the function ■0^1 is 
compactly supported and smooth in the spatial variables, its Fourier transform 


has rapid decay. This shows that the function (3.20) is indeed square integrable. 


Using Plancherel, we see that condition (a) in Definition 2.3 is satisfied. Moreover, 
the operator T is simply the operator of multiplication by v^, so that condition (b) 
obviously holds. This again shows the weak mass oscillation property. 
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In order to prove the strong mass oscillation property, we need to compute the inner 
product <p'0|p(/)>. To this end, we first write this inner product in momentum space 
as 


<p'0|p^> = 


d^k 

d‘^k 


4vr^ XI {m) 


1 




4m^ 


010 








Reparametrizing the fe^^-integral as an integral over m, we obtain 


^7Vm(fe) I (Hm) 7Vm(fe)71• (3-21) 


Estimating the inner product with the Schwarz inequality and applying Plancherel’s 
theorem, one finds 


|<pV’|p(?i>| < J^dm 


i(^)ll d k 27r J ||'0m||m ||0m||m dfU . 


Thus the inequality (2.15) holds. 


□ 


4. The Mass Oscillation Property in Minkowski Space with External 

Potential 

4.1. Proof of the Weak Mass Oscillation Property. In this section, we prove 
the following theorem. 

Theorem 4.1. Assume that the time-dependent external potential 23 is smooth and 


decays faster than quadratically for large times in the sense that (1.7) holds for suitable 


constants c,e > 0. Then the Dirac operator D = i<^ -\-‘B has the weak mass oscillation 
property. 

We expect that this theorem could be improved by weakening the decay assumptions 
on the potential. However, this would require refinements of our methods which would 
go beyond the scope of this paper. Also, using that Dirac solutions dissipate, the 
pointwise decay in time could probably be replaced or partially compensated by suit¬ 
able spatial decay assumptions. Moreover, one could probably refine the result of the 
above theorem by working with other norms (like weighted C^- or Sobolev norms). 

The main step is the following basic estimate, which is the analog of Lemma |3.1| in 
the presence of an external potential. 


Proposition 4.2. Under the decay assumptions (1.7) on the external potential 23, 


there are constants c, e > 0 such that for every family G TC“ of solutions of the 


Dirac equation (2.3) with varying mass, 


|(pV')ML < . , sup ||(5^iAm)|t=o 

r “T Ul mdl ■ - 


=u 1114/2.2 


(4.1) 


b=0 


We first show that this proposition implies the weak mass oscillation property. 

Proof of Theorem \4.1\ un der the assumption that Proposition \4.S\ holds. In order to de¬ 
rive the inequality (2.13), we begin with the estimate 


\<p'ip\p(j)>\ < — 


(pV’lt I P</>lt)lt dt < sup||p(/)|t| 

ieIR 


Ipi^Ulldt. 
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The last integral is finite by Proposition |4.2[ The supremum can be bounded by the 
Hilbert space norm using the Holder inequality, 


MltWt = 


^m\t dm 


< 


\4>m\t\\^ dm < ( 


\Ji 


bm\t\\t dm 


= \n 


giving ( |2 13D. 

Using (2.2), the Dirac operator + 23 is formally self-adjoint with respect to the 
inner product <.|.>. Therefore, the identity (2.14) can be obtained just as in (3.6) 
by integrating the Dirac operator in space-time by parts, noting that we do not get 
boundary terms in view of the time decay in Proposition 4.2 □ 


The remainder of this section is devoted to the proof of Proposition 4.2 One 


ingredient is the Lippmann-Schwinger equation (2.20), 

rt 


V'mit = V'mk=0 + i / (7°® 1 pm)\^ dr . 


(4.2) 


'0 


Since the first summand of this equation is controlled by Lemma 3.1, it remains to 
estimate the second summand. Again using (3.15) and integrating by parts with 
respect to the mass, we obtain 

1 


'lpm)\^dm = 


{t - t)2 


and thus 




< 


(7°® V’m) dm 

c\I\ 


it - r)2 


sup 

VF2'- 


m&I 


b=0 


< 


sup 

U mei ^^0 


\T\\^r2,2 ■ 


We now bound 23 (r) with the help of 0 and estimate the Sobolev norm ||9^V'm|T||y[/2.2 
at time r by means of Lemma |A.1| proved in Appendix This gives rise to the in¬ 
equality 

2 


'lpm)\^dm 


c^C\I\ l + |r |2 


t I + |r|2+£ 


sup 


'mU=0||^y2,2 ) 


b=0 


which yields the desired decay provided that r and t are not close to each other. More 
precisely, we shall apply this inequality in the case |r| < |t|/2. Then the estimate 
simplifies to 


dm 
C 


< 


,2n I \\9ti^m\t=o\\w2.2 if kl < \t\/2 

I {1 -\- \T\ ) rnGl 


(4.3) 


with a new constant C > 0. In the remaining case |t| > |t|/2, we use the unitarity 


of to obtain 


(7°® V'm) It dm 


< \I\ l®('^)lco sup I 

t m£l 
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Applying 0 together with the inequality |r| > |t|/ 2 , this gives 

(7°® V'm) |r dm <-^ 


t2+e 


sup I 

mGl 


if |r| > |t |/2 . 


(4.4) 


This again decays for large t because r is close to t and decays for large r. 

Comparing (4.3) and (4.4), we find that the inequality in (|4.3l) even holds for all r. 


Thus integrating this inequality over r G [0,t], we obtain the following estimate for 
the second summand in (4.2), 


I dm J ( 7 °® i/jm) \t dr 


< 


(jf 

^sup '^\\di,ip\t=o\\^2,2 


(where C' > 0 is a new constant). Combining this inequality with the estimate (3.2) 
of the first summand in (|4.2[ ), we obtain the desired inequality (4.1). This concludes 
the proof of Proposition |4.2[ 

4.2. Proof of the Strong Mass Oscillation Property. In this section, we prove 
the following result. 

Theorem 4.3. Assume that the weak mass oscillation property holds and that the 
external potential 13 satisfies the condition 

/ OO 

\'B{t)\co dr < 00 . (4.5) 


Then the Dirac operator D = it^ + iB has the strong mass oscillation property. 


Combining this theorem with Theorem 4.1 one immediately obtains Theorem 1.2 


For the proof we shall derive an explicit formula for the fermionic signature op¬ 
erator (Proposition 4.4). This formula is obtained by comparing the dynamics in 
the presence of the external potential with that in the Minkowski vacuum using the 
Lippmann-Schwinger equation, and by employing distributional relations for products 


of fundamental solutions and Green’s functions (Lemma 4.7). 

We first return to the formula (3.21) in the Minkowski vacuum. Applying Plancherel’s 


theorem and using (2.4), we conclude that 


where 


<pfi\p(j)> = I $m{k) dm , 


Q l^ + m 0 kfi + m 0 

^m[k) := 2 ^ 7 '^ = - 7 '' 


/c®=±aj(Ai) 


2uj{k) 


uj{k) 


(4.6) 

(4.7) 


Comparing (4.6) with (2.16), one sees that the matrix §rn{k) is indeed the fermionic 


signature operator, considered as a multiplication operator in momentum space. By 
direct computation, one verifies that the matrix §rn{k) has eigenvalues ± 1 . 

In order to compare the dynamics in the presence of the external potential with that 
in the Minkowski vacuum, we work with the Hamiltonian formulation. We decompose 


the Dirac Hamiltonian (2.19) into the Hamiltonian in the Minkowski vacuum (2.21) 
plus a potential, 

H = H + V with V := - 7 °® . 
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Proposition 4.4. Assume that the potential 23 satisfies the condition (4.5). Then for 
every ^ !H°°, 

Kpfilpfiy = / {firn I Sm 4>m)m dm , (4.8) 


where §>ra '■ 'dim 
at time to by 


dim are bounded linear operators which act on the wave functions 


Sm — Sm, „ 


1 

+ 2 


s(t - to) [s™ V(t) V(t) Sm dt (4.9) 


—oo 
oo roo 



+ 


to Jto 


do rto 


’ —oo J —oo 


V(t} Sm V{t') dt dt' (4.10) 


(and Sm is again the fermionic signature operator of the vacuum (4.7)). 


Before entering the proof of this proposition, it is instructive to verify that the above 
formula for Sm does not depend on the choice of to. 


Remark 4.5. (Independence of Sm on to) Our strategy is to differentiate the 
above formula for Sm with respect to to and to verify that we obtain zero. We first 
observe that taking a solution fim £ dim of the Dirac equation in the presence of 23, 
evaluating at time to and applying the time evolution operator Um° gives fim at time t, 
i.e. Um''(t)m\to = fim\t- Differentiating with respect to to yields 

^to Um ° 4’m I to = 0 . 


The situation is different when one considers the time evolution operator of the vac¬ 
uum. Namely, in the expression Um° fimUo^ the wave function fim satisfies the Dirac 
equation (idt — H)(j)m = whereas the time evolution operator solves the Dirac 

equation with V = 0. As a consequence, 

dtoUl^^^m\to = -iU(d''{Vfim)\to ■ 

Using these formulas together with l straightforward computation 

gives 


dto {"fim I ( |4^ fim) I 


^0 




^0 


- i{fim I [SmiV] cj)m)\to 

- ^ (-2) (V’m I (Sm ^(to) - "i^(to) Sm) 4>m) 1^^ 

- A / e(t - to) ((-iV(to)) fim I Sm V(t) U(do (j)^')\^^dt 
i 

+ A / <t- to) ifim I uy V{t) Sm {-iV{to)) fim) 

1 

- A / <t' - to) {fim I V(to) Sm V{t') 1]^^^ dt' 

I 

- A / - to) {fi’m I f/Do y7(iQ) (j)^') dt , 

^ J —OO 


where for notational simplicity we here omitted the restrictions It,, for the solutions fim 
and (fm- Adding the terms gives zero. 0 


The remainder of this section is devoted to the proof of Proposition |4.4[ Our strategy 
is to combine the Lippmann-Schwinger equation with estimates in momentum space. 
We begin with two technical lemmas. 
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Lemma 4.6. Assume that the external potential "B satisfies condition (4.5). For 
any to £ we denote the characteristic functions in the future respectively past of 
this hypersurface t = t^ by (i.e. = 0(±(x° — to)), where 0 is the Heavi¬ 
side function). Then for any the wave function 

is a well-defined vector in and 

1 


2-k 


Xw{t) I'B(t)|co dr . 


Proof. Using the integral kernel representation (2.8) and (2.9) together with the fact 
that the time evolution in the vacuum is unitary, we obtain 


27r 


km{{to,.), (r,y)) {xtgBfirn){T,y) d^y 


to 


= = h\xtMm)\r\l < |®(r)|cO 


Integrating over r and using (4.5) gives the result. 

The following lemma is proved in |18l Eqs. (2.13)-(2.17)]. 


□ 


Lemma 4.7. In the Minkowski vacuum, the fundamental solution km and the Green’s 
function Sm defined by 

1 


Sm. - 1 + s)),) 


(4.11) 


satisfy the distributional relations in the mass parameters m and m' 
km km' = 5{m- m') pm 


km Sm' — Sm' km — 


PP 


m — m 


PP 


^m. Stri/ — 


7 (Sm - Sm') + <5(m - m') Pm , 


m — m' 

where PP denotes the principal part, and pm is the distribution 

Pm{k) = (1^ + m) 6{k‘^ - w?) . 


(4.12) 


Proof of Proposition 4-4 ■ Let if G be a family of solutions of the Dirac equation 
for varying mass. We denote the boundary values at time to by i/’m := tfm\to- Then 
we can write the Lippmann-Schwinger equation (2.20) as 

ifmlt = + i [ '>Pm) dr . 

Jto 

We now bring this equation into a more useful form. Expressing the time evolution 
operator with the help of (2.9) in terms of the fundamental solution, we obtain 

ifm{x) = 2 tt km{x, (^0,2/)) 7°V’m(io,27) d^y 


+ 27ri dy^ d^y km{x,y){B fimfiv) ■ 


'to 


Applying (2.5) and using that the advanced and retarded Green’s functions are sup¬ 
ported in the future and past light cones, respectively, we can rewrite the last integral 
in terms of the advanced and retarded Green’s functions, 

ifm = 2-Kkm{'l^dt^,lpm) - S^{xto'^'fi'm) “ , 
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where 5to{x) := 5{to — x^) is the Dirac distribution supported on the hypersurface x^ = 
to- Next, we express the advanced and retarded Green’s functions in terms of the 
Green’s function (4.11): According to (2.5), we have the relations 

Sm — iirkm — + iirkm 

and thus 

= kmOm - with - •- ./.O 


gm ■= 271 + ill eto'B'iprn , (4.13) 

where tto is the step function 

eto(x) := e(x° - to) 

(and we omitted the brackets in expressions like kmOm = km{gm))- Note that the 
expression kmgm is well-defined according to Lemma 4.6 We also remark that by 
applying the operator — m) to the distribution gm in d^ , one immediately 
verifies that indeed satisfies the Dirac equation (z^ — m)V’m = — 


Now we can compute the inner product with the help of Lemma 4.7 


Namely, using (4.13), 


<p'i/’|p'0> — II ^kmgm | km'gm' dm dui 


Ixl 


— i^^gm I Pmgm^ I dm 

^ I km'gm'^ ^kmgm \ 

-7 <‘Biljm I (sm - 'Sm')®V’m'>) dm dm' 


Ixl m-m' 


Note that this computation is mathematically well-defined in the distributional sense 
because z/'m and gm are smooth and compactly supported in the mass parameter m. 
Employing the explicit formula for gm in (4.13), we obtain 

^pz/^lpz/^> = i^^gm I Pmgm^ 4“ H \ dm . 


Comparing 


with (4.12) and taking into account that the operator §m defined 


by (4.7) gives a minus sign for the states of negative frequency, we get 

Pm — Sm km ■ 


Using this identity together with Proposition |2.2| in the vacuum yields the relations 

^9m I Pm9m^ — {km9m \ Sm ^m5m)|io 

^^V’m I Pm‘^'4^m,^ — (^m'^V’m | ^m fcm®V’m)|to • 

We finally apply Proposition |2.1| to obtain the representation 


<pz/^|p'i/’> — ^(^m I Sm /i-m)|to 4~ ^ (^m^V’m | Sm ^m®V’m)|zo^ dm , (4.14) 


where 


km ■ — V’m 4“ in kmif-tQ^'^^m') ■ 

Comparing (4.8) with ( 4.14[ ), we get 

(V’m I Sm '^m)m — {hm I Sm km)\tQ 4" VT (/Cm®V’m | Sm ^m'®V’m)|to 
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Expressing the operators km according to (2.9) by the time evolution operator and 
writing ij^m in terms of the initial data as 


^m\t = 


we obtain 


{'Ipm I Sm 'lpm)m = ('0|SmV’)|to “ X / ~ I Ln 


1 


1 

+ T 


it-to) U^’^°^|J\§mi)\tJt 

e{t - to) e{t' - to) ([/*“’* V(t) i I V{t') V’) I* dt dt' 


^t,to ^ I 1^ ^ 


Rearranging the terms and polarizing gives the result. 


□ 


Proof of Theorem \4.3[ Since the time evolution operators are unitary and the opera¬ 
tors §m have norm one (see (4.7)), the representation (4.9) and (4.10) gives rise to the 


following estimate for the sup-norm of Sm, 
||Sm|| < 1 + [ \'^it)\c°dt + 


\m\co \W)\co dtdt' . 


The decay assumption (4.5) implies that the sup-norm of S', 
in m 


is bounded uniformly 


Using this fact in (4.8) gives the inequality (2.15), thereby establishing the 

□ 


strong mass oscillation property. 


We finally remark that the uniqueness statement in Theorem 2.5 implies that (4.9) 


and (4.10) yields an explicit representation of the fermionic signature operator in the 


presence of a time-dependent external potential. 


5. Hadamard Form of the Fermionic Projector 


In this section, we will prove Theorem 1.3 


In preparation, we derive so-called fre¬ 
quency splitting estimates which give control of the “mixing” of the positive and nega¬ 
tive frequencies in the solutions of the Dirac equation as caused by the time-dependent 
external potential (Theorem 5.1). Based on these estimates, we will complete the proof 
of Theorem O at the end of Section 15.21 


5.1. Frequency Mixing Estimates. For the following constructions, we again choose 
the hypersurface jV := at some given time to- Moreover, we always fix the mass 
parameter m > 0. Since we are no longer considering families of solutions, for ease in 
notation we omit the index m at the Dirac wave functions, the scalar products and 
the corresponding norms. We also identify the solution space “Km with the Hilbert 
space Jitg of square integrable wave functions on W. On TC^p, we can act with the 
Hamiltonian H of the vacuum, and using the above identihcation, the operator H 
becomes an operator on Jim (which clearly depends on the choice of to). 

We work with a so-called frequency splitting with respect to the vacuum dynamics. 
To this end, we decompose the Hilbert space Jim as 

Jim = Jim® Jim with Ji^ = , 
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where a-re the characteristic functions 

:= X[0,oo) and x" := X(-oo,o) • (5-1) 

For convenience, we write this decomposition in components and use a block matrix 
notation for operators, i.e. 

^ and d = jz'j , 

where A^, = x^iH)Ax^'(H) and s,s' G {±}. 

The representation in Proposition |4 . 4| makes it possible to let the fermionic signature 
operator act on the Hilbert space "Km (for fixed m). We decompose this operator 
with respect to the above frequency splitting, 

Sm = + AS , where S° := S^ + Si and AS := Sl + Si,! . 

Thus the operator S° maps positive to positive and negative to negative frequencies. 
The operator AS, on the other hand, mixes positive and negative frequencies. In 
the next theorem, it is shown under a suitable smallness assumption on "B that the 
operators x^(Sm) coincide with the projections x^{H), up to smooth contributions. 
The main task in the proof is to control the “frequency mixing” as described by the 
operator AS. 

Theorem 5.1. Under the assumptions of Theorem 
the representations 

(Sm) = X^{H) + ^.(L [Im - A)-1 AS (S° - A)-^ dA , (5.2) 

JqBx (± 1 ) 

2 

where the contour integral is an integral operator with a smooth integral kernel. 


1.3. the operators x^(Sm) have 


Here Bi denotes the open ball of radius 1/2. The operator (Sm — A) ^ is also referred 
2 

to as the resolvent of §m- 

This theorem will be proved in several steps. We begin with a preparatory lemma. 


Lemma 5.2. Under the assumptions (1.7) and (1.14), the spectrum o/S° is located 
in the set 


ct(S°) C 


3 

2 ’ 


U 


1 3 
2 ’ 2 


Moreover, 


X^(S°)=X^(i^), 


±/ 


(5.3) 

(5.4) 


and the operators x^(Sm) have the representations (5.2). 


Proof. Since the subspaces are invariant under the action of S°, our task is to 
show that the spectrum of S°|j{;± is positive and negative, respectively. This state¬ 
ment would certainly be true if we replaced S° by Sm, because the operator Sm has 
the eigenvalues ±1 with as the corresponding eigenspaces. Estimating the repre¬ 
sentation in Proposition 4.4 with the Schwarz inequality, we obtain 

r.2' 


< {c+— ] \\fj\\\\4>\\ with c:= |®(r)|codr 
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Using the assumption (1.14), we conclude that 

J Ill/'ll 1101 


for all '4),^ ^ CKf, 


Standard estimates on the continuity of the spectrum (see for example [33 §IV.3]) 
yield that the spectrum of differs by that of the operator Sm at most by 1/2. This 
gives (5.3) and (5.4). 

In order to prove the representation ( |5.2| ), we take the resolvent identity 
(Sm - A)-1 = (S° - A)-1 - (Sm - A)-1 AS (S° - A)-1 , 


form the contour integral and apply (5.4). This gives the result. 


□ 


The next lemma relates the smoothness of an integral kernel to the boundedness of 
the product of the operator with powers of the vacuum Hamiltonian. 

Lemma 5.3. Let A G L(lKm) be an operator which maps smooth functions to smooth 
functions and has the property that for all p, g G N, the operator product 


H^ARP : C^(^,SM) ^ C’^(J(,SM) 


(5.5) 


extends to a bounded linear operator on Jim- Then, considering A as an operator 
on TCm, this operator can he represented as an integral operator with a smooth integral 
kernel, i.e. 

(A4j)(x) = f A(x, (to,y))'ijj(tQ,y) d^y with AgC°^(MxM). 

Proof. Since in momentum space, the square of the Hamiltonian takes the form 
H(k)^ = ^ 7 ^( 7 ^ + = ( — 7 ^ + m) ( 7 ^ + m) = \kf + , 

the wave function 0 defined by 


i>(k) := 


1 


12 _|_ 


„ikxo 


for a constant spinor H and xq G satisfies the equation 

'ij:(x) = 6^(x — xq) H . 

Moreover, one verifies immediately that 0 G “Kta is square-integrable. Using the last 
equation together with (|5.5|), we conclude that 


H^A(5^(x - aTo) ") = H^AH^fj G Rto • 

Since q is arbitrary, it follows that A has an integral representation in the spatial 
variables, 

(A4>)(x) = f A(x, y) 7 *^ 0(y) d^y with A G C°°(J( x A) . 

Jm 

We now extend this integral kernel io Ji x Ji hy solving the Cauchy problem in the 
variables x and y. This preserves smoothness by the global existence and regularity 
results for linear hyperbolic equations, giving the result. □ 







28 


F. FINSTER, S. MURRO AND C. ROKEN 


Lemma 5.4. Under the assumptions of Theorem for all p € N the iterated com¬ 
mutator 

gip) 


H, 


p factors 


is a bounded operator on TCm- 


Proof. In the vacuum, the Hamiltonian clearly commutes with the time evolution 
operator, 

[H,[/^*']=0. (5.6) 

In order to derive a corresponding commutator relation in the presence of the external 
potential, one must take into account that H is time-dependent. For ease in notation, 
we do not write out this dependence, but instead understand that the Hamiltonian is 
to be evaluated at the correct time, i.e. 


H = H{t') and HU^/ = H{t)U^/ 


Then 


{idt - H) {H Utf - Ut^'H) = iH and 


H Utit' - Utit'H 


t=t' 


= 0 


(here and in what follows the dot denotes the partial derivative with respect to t). 
Solving the corresponding Cauchy problem gives 

[H,U^J] = dr. (5.7) 


In order to compute the commutator of H with the operator products in (4.9) 


and (4.10), we first differentiate the expression Um'^VUm with respect to t 
idt i V (7*’^') = V f/*’*' -b V H 17*’*' - 7/*"’* H V t/*’*' 


jtf' 


(5.8) 


Moreover, using the commutation relations (5.6) and (5.7), we obtain 
= t/r H V 7/4*' - 77^’* V H iJl^' + 77^’* V [#, 77^*'] 

= i77^'’* V 777*' _ JjtUt y jjrJ ^ 

Jt' 

where in the last step we applied ( |5.8[ ). It follows that 

[77, 77^'’* V ^^*'] = 77 (77^'’*W7*') - (77^'’*W^*') H + (77^'’*V^4*') V 

= 777^'’* V 777*' + (77^'7v777*') V - idt V 777*') + f' jjtPt ^ fjt^r p jjrj 

Jt' 

Proceeding in this way, one can calculate the commutator of 77 with all the terms 
in (4.9) and (4.10). We write the result symbolically as 

[77,S^] = S«, 

where 8^*^ is a bounded operator. Higher commutators can be computed inductively, 
giving the result. □ 
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We point out that this lemma only makes a statement on the iterative commutators. 
Expressions like [H^, S^] or §rn will not be bounded operators in general. How¬ 
ever, the next lemma shows that the operator AS has the remarkable property that 
multiplying by powers of H from the left and/or right again gives a bounded operator. 


Lemma 5.5. Under the assumptions of Theorem 
H‘> AS is a bounded operator on ■ 


jor oh p, g G Jn U the product 


Proof. We only consider the products S_,_ HP because the operator Sl can be treated 
similarly. Multiplying (5.7) from the left and right by the resolvent of H, we obtain 

[{H - = -{H - /r)-i S« {H - /r)-i . 

Writing the result of Lemma |5.4| as 

[/7,S(P)] = §(P+i) with S(P+^) G L(Jf) 

yields more generally the commutation relations 

[{H - fi)-\§^P^] =-{H - fi)-^§^P+^'> {H - forpGN. (5.9) 

Choosing a contour 7 which encloses the interval (— 00 , —m] as shown in Figure 
one finds 

= dfi 

= §Hx-{H) x+(H) + ^ ^ /r (77 - /r)-i - /r)-^ x+(i7) d^^ 

= ^ ^ h (^ - h)-' - p)-^ y+(H) dp , 

where in the last step we used that = 0. In order to show that this 

operator product is bounded, it is useful to employ the spectral theorem for 77, which 
we write as 


f{H) = / 

^R\[—m,m] 

where dEx is the spectral measure of 77. This gives 

1 


77 SI = 


T 


(± 

\ 27ri J.y X - p X' - p 
X 


/(A) dEx , 

ves 

X+(A') dEx ) S(^) dEy dp 


(5.10) 


A-A' 


X"(A) x+(A')dEAS(^) dEy. 


(5.11) 


Note that the term A — A' is bounded away from zero. Thus the factor A/(A — A') is 
bounded, showing that the operator H§f_ is in L(TCm). 
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This method can be iterated. To this end, we first rewrite the product with com¬ 
mutators, 

H‘1 s; = x" (H) {H- X- {H)Y s™ x^{H) 

= X- (H) [h-,[H-,...,[H-,§]---]]x^ (H) , 

where we used the abbreviation H~ := Hx~{H). Multiplying from the right by 
we can commute factors H'^ := H x'^{H) to the left to obtain 

§,-HP = i-ir x-iH) ■■■]]■■■] X^iH) . 

^ ■v* 

p factors q factors 

Representing each factor by a contour integral, one can compute the commutators 
inductively with the help (5.9). Applying the spectral theorem (5.10) to the left and 
right of the resulting factor yields a constant times the expression 


X-(A) x^{^')dExS^P+<i'> dEy 

/ii d^i 


71 




fJ'P+q dfip+q 


7p+9 


(A d-p+q){^' d-p+q) 


Carrying out the contour integrals with residues, we obtain similar to (5.11) an ex¬ 
pression of the form 




fix, A') x"(A) x+(A') dEx dEx' 

with a bounded function /. This concludes the proof. 


□ 


Proof of Theorem 5.1. It remains to be shown that the contour integral in (5.2) has 
a smooth integral kernel. To this end, we multiply the integrand from the left by 
and from the right by and commute the factors H iteratively to the inside. More 
precisely, we use the formula 

- A)-‘ = V [if,..., [if, {S„ - A)-‘]. ■ ■ ] if’"“ 

(note that the sum is telescopic; here we use the convention that the summand for a = 0 
is simply (S^ — X)~^H^). Hence 

H'' (§™ - A)-i AS (S° - A)-i HP 

q p 


H,...,[H,(S™-A)-i]--- 


a=0 6=0 


AS HP 


-b 


-X)-\H],...,H 


a factors 


b factors 


According to Lemma 5.5, the intermediate product H^ AS HP " is a bounded op 


erator. Moreover, the commutators can be computed inductively with the help of 
Lemma 15.41 and the formula 

[H, (Sm - A-1)] = -(S™ - A-l) [H,~Sm] i§m " A'^) 

(and similarly for S^). This gives operators which are all bounded for A G cli?i(±l). 
Since the integration contour is compact, the result follows. □ 
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5.2. Proof of the Hadamard Form. Relying on the frequency mixing estimates 
of the previous section, we can now give the proof of Theorem 1.3 Recall that the 


fermionic projector is given by (see (1.9)) 


P = -X (Sm) K 


(5.12) 


where we again used the short notation (113. Here again the operator x (Sm) acts on 
the solution space TCm of the Dirac equation, which can be identified with the space 
of square integrable wave functions at time to (see the beginning of Section 5.1). For 


the following arguments, it is important to note that this identification can be made 
at any time to¬ 
la order to prove that the bi-distribution corresponding to P is of Hadamard form, 
we compare the fermionic projectors for three different Dirac operators and use the 
theorem on the propagation of singularities in [38]. More precisely, we consider the 
following three fermionic projectors; 

(1) The fermionic projector in the Minkowski vacuum. 

(2) The fermionic projector P in the presence of the external potential 

B(x) :=r?(x°)®(x), 

where r/ > 0 is a smooth function with ??|(_oo,0) = 0 ??|(i,oo) = 1- 

(3) The fermionic projector P in the presence of the external potential 13(x). 

The potential ‘B vanishes for negative times, whereas for times > 1 it coincides 
with “B. Thus it smoothly interpolates between the dynamics with and without external 
potential. The specific form of the potential B in the transition region 0 < x^ < 1 is 
of no relevance for our arguments. 


In the Minkowski vacuum, the relation (5.12) gives the usual two-point function 
composed of all negative-frequency solutions of the Dirac equation. It is therefore 
obvious that the bi-distribution P™'^(x,y) is of Hadamard form. 

We now compare with P. To this end, we choose an arbitrary time to < 0. 


Then, applying the result of Theorem 5.1 to (5.12), we get 


pvac ^ p ^ p (smooth) 


where km is the causal fundamental solution in the presence of the potential B. Since B 
vanishes in a neighborhood of the Cauchy surface at time to, we conclude that P™'^ 
and P coincide in this neighborhood up to a smooth contribution. It follows that 
also P(x, y) is of Hadamard form in this neighborhood. Using the theorem on the 
propagation of singularities |38[ Theorem 5.5], we conclude that P{x, y) is of Hadamard 
form for all x, y G . 

Next, we compare P with P. Thus we choose an arbitrary time to > 1. Using again 
the result of Theorem 5.1 in (5.12), we obtain 

P = — (H) km + (smooth) and 


P = —X (H) km + (smooth) 


(where the smooth contributions may of course be different). Since B and B coincide in 
a neighborhood of the Cauchy surface at time to , we infer that P and P coincide in this 
neighborhood up to a smooth contribution. As a consequence, P{x, y) is of Hadamard 
form in this neighborhood. Again applying |38l Theorem 5.5], it follows that P{x,y) 
is of Hadamard form for all x, y G .y^. This concludes the proof of Theorem |1.3[ 
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6. Quantum Fields and the Hadamard State 


In preparation of defining the CAR algebra, we denote the co-spinor bundle by S*Ji 
(thus the fiber S^M is the dual space of SxM). On the smooth and compactly sup¬ 
ported sections of S*Ji, we introduce the dual of the Dirac operator D* by 

V* : , 

[ {(V*g){f)){x)d^x= [ {g{Vf)){x)d^x for all / E Co~(^, 5^) . 

JjU, Jjt 

Moreover, we define the space of pairs of spinorial test functions by 

with the topology induced by the family of seminorms 

l(/, 5 )lc'= := sup |5^/(x)| sup \d’'g{y)\ 

xdJi y&jii 

(and |.| is any norm on the spinors and co-spinors, respectively). Next, we define the 
anti-linear involution map 

FiD^D, r{f®g):=g*®r, (6.1) 

where the star again denotes the adjoint with respect to the spin scalar product, i.e. 


* : C^{Jl,S*M)^ C^{Jl,SJl), ^g*\(py = g{(P) 


(we remark that the adjoint spinor f * can be identified with the adjoint spinor usually 
written as / = /^ 7 *^). Finally, we introduce an inner product on D, 


(.|.)s) : 2) X DC , {f ® g \ aeib )^ = < f \ krna > + < b *\ kmg *> , (6.2) 

where <.|.> is again the space-time inner product 
with our overall convention that brackets like (.|.), (.|.)m and <.|.> are sesquilinear 
forms, meaning that the first argument always involves complex conjugation). Apply¬ 
ing Proposition |2.2[ we can write this inner product as 


( 2 . 1 ) (this notation is consistent 


(/ © 5 I a © 6 ) = (kmf I kma)m + (kmb* I km g*)m , 


showing that it is indeed positive definite. Forming the completion, we thus obtain the 
Hilbert space (Rs > (-I-)®)- A short computation reveals that the involution F preserves 
the scalar product (.j.)® and is therefore anti-unitary, i.e. 


F^ = 1 and {Th \ Th)^ = {h\h)^ for all h,h € fhx) • (6-3) 

Definition 6.1. The field algebra T' is the unital *-algebra generated by the abstract 
element B{h) with h £ T), together with the following relations for all f ® g £ “T) 
and a, (3 £ C: 

(i) Linearity: B{af ® g + m® jdn) = aB{f © 5 ) + f3B{m © n) 

(ii) Hermiticity: B(f © g)* = B{r{f © g)) 

(hi) Dynamics: B{^(T) — m)f © {V* — m)g) = 0 for all f ® g £ T) 

(iv) Canonical anti-commutation relations (CARs) 

{B{f © gT,B{m ® n)] = {f ® g\m ® n)^ ■ Ijr . 
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We define smeared field operators by 

:=S(0©5) 

Then these held operators satisfy 


and T*(/) :=i?(/©0), 


{T(5), '&*(/)} = {B{0 © g), Bif © 0)} {S(0 © g),B{T{f © 0))*} 


(r(/ © 0) I 0 © 5)2, ^ (o © /* I 0 © 5)2, ^ <g* \kmf>, 

giving rise to the usual anti-commutation relations ( 1.17[ ). 

The next step is to show that the fermionic projector induces a pure, quasi-free 
state on the held algebra. Let us recall a few basics. By dehnition, a state w is a linear 
functional on B which is positive and normalized, i.e. 


u{B*B) > 0 for all B G T" and = 1 ■ 

By linearity, it suffices to specify u on the monomials. This gives rise to the n-point 
distributions a;„ G (D"')' dehned by 

^n(hli • • • 1 hn) .— ^{^B(h\^ • • • B(^hn')'j , 

also referred to as the n-point functions. A state induces the usual Fock representation 
of the held algebra: 


Theorem 6.2. (GNS construction) Let uj be a state on a unital *-algebra T. Then 
there exists a Hilbert space {Hpock, (■, ■)Fock) o-nd a dense subspace W C Jipock as well as 
a representation vr : T” —)• L(>V) and a unit vector n G W such that oj = (fl, 7r(.)0) po^k 
and W = 7r{iF)Q. The GNS triple (>V,7r,n) is determined up to unitary equivalence. 


Among all possible states, a distinguished role is played by the quasi-free states, for 
which the n-point functions are all determined by the two-point functions: 


Definition 6.3. A state cu : T" —)• C is called quasi-free if the n-point functions vanish 
for odd n, while for even n, one has 

nl2 

UJn{hl,...,hn) = J]] n ^2(/ia{2i-l),^a{2*)) , 

i=l 

where 5^ denotes the set of ordered permutations of n elements, i.e. 

Tl 

a{2i — 1) < a{2i) for all i = 1,..., — 

Tl — 2 

a{2i — 1) < a{2i + 1) for all i = 1,..., —-— . 

For the construction of a quasi-free state from the fermionic projector, we rely on the 
following result due to Araki [U Lemmas 3.2 and 3.3.]: 

Lemma 6.4. Let R be a bounded symmetric operator on (IKs), (.|.)j)) with the following 
properties 

(a) R + Fi2F = 1, 

(b) 0 < i? = ii* < 1 . 

Then there exists a unique quasi-free state uj on H such that 

uj{B{h)*B{h)) = {h\Rh)s for all h,h e . 


(6.4) 
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Moreover, Araki proves that if ii is a projection operator, then the state oj is pure 
(see [U Lemma 4.3]). 

Thus our task is to construct an operator R with the above properties using the 
fermionic signature operator. As will become clear below, the correct choice is to 
define R implicitly by 

[f ® 9\R{a®b))^ = (fcm/l X"(Sm) + {k^b* I kraQ*)^ , (6.5) 

giving a bounded and symmetric operator. Moreover, since x^(Sm) are projection 
operators, we know that 


0 < (/ © £/ I i? (/ © £/)) < {kmf I kmf)^ + {kmg* \ kruG*)^ = (/ © 5 I / © ff)2) , 

showing that the condition (b) in Lemma |6.4| holds. Next, using that L is an anti¬ 
unitary involution ( |6.3[ ) and R is symmetric, it follows that 

(/©5|ri?r(a©6))^ = (i?r(a©6)|r(/©9))^ = (r(a©6)|i?r(/©(7))^. 

Now we can apply and ( |6.5[ ) to obtain 

(/ © <71 TRY (a © 6)) = {b* (B a* \ R {g* © f))^ 

= (kmb* I X“(Sm) kmg*)^ + (kmf \ X’^(Sm) • (6.6) 


Adding ( |6.5[ ) and ( |6.6[ ), we get 

{f®g\{R + YRT) (a © b))^ = (kmf \ kmo)^ + [kmb* \ kmg*)^ = [f ®g\a®b)^, 


proving that the condition (a) in Lemma |6.4| is satisfied. Thus Lemma |6.4| applies, 
giving a quasi-free state iv with the property (6.4). Moreover, as is verified by direct 
computation, i? is a projection operator. Therefore, the state u is pure. 

We finally calculate the two-point function. Beginning with the computation 


u;{^{g)^*{f)) =uj[B[0(Bg)B{f(B0)) u;(^B{T{ 0 (B g))* B{f (B 0)) 

® oj{B{g* © 0))* B{f © 0)) ® (5* © 0 I i? (/ © 0))^ 
® (km9* I X~(§m) kmf)^ ^ <9* \ X~ i^m) kmf> 


(where in (*) we used the property (ii) in Definition 6.1), we can apply the definition 
of the fermionic projector (2.17) to obtain 


a;(T( 5 )T*(/)) = -<g*\Pf> = - jj g{x)P[x,y)f{y)(fx(fy. 


This concludes the proof of Theorem 1.4 


Appendix A. Uniform L^-Estimates of Derivatives of Dirac Solutions 
We now derive a few estimates which will be needed for the proof of the mass 


oscillation property in Section 4.2 We use standard methods of the theory of partial 


differential equations and adapt them to the Dirac equation. In generalization of (3.3), 
we denote the spatial Sobolev norms by 




E 

a with \a\ < a 


|V“())(f)r 
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Lemma A.l. We are given two non-negative integers a and b as well as a smooth 
time-dependent potential 23. In the case a > 0 and b > 0, we assume furthermore that 
the spatial derivatives of 23 decay faster than linearly for large times in the sense that 

for suitable constants c, e > 0. Then there is a constant C = C{c, e, a, b) such that every 
family of solutions ^|J G !K°° of the Dirac equation (2.3) for varying mass parameter 
can be estimated for all times in terms of the boundary values at t = 0 by 

b 

\\dL'^m\t\\wa ,2 <C{1 + |t|^) \\d^'lpm\t=o\\wa ,2 ■ 

p=0 

Proof. We choose a multi-index a of length a := |a| and a non-negative integer b. 
Differentiating the Dirac equation (2.3) with respect to the mass parameter and to the 
spatial variables gives 

+ m) = bS7‘^^t^^Pm - V“(iB ® • 

Introducing the abbreviations 

H := and := 6 Vm - V“(23 + 23 , 

we rewrite this equation as the inhomogeneous Dirac equation 

{V — m)E = <f>. 

A calculation similar to current conservation yields 

—idj^E\^^E>~ = ^{V — m)E \ -<H | {V — m)E>- = -<(/)|H>- <E\(I)'^ . 

Integrating over the equal time hypersurfaces and using the Schwarz inequality, we 
obtain 




and thus 


d. 


< mt\ 


Substituting the specific forms of E and </> and using the Schwarz and triangle inequal¬ 
ities, we obtain the estimate 


9d|V“5iV'™|dL <b\\V^dl^f;m\t\i+ca\V'B{t)\ca-i 


jaab—l 


aab 


liy“ 


(A.2) 


where we used the notation (1.8). 

We now proceed inductively in the maximal total order a -|- 6 of the derivatives. 
In the case a = b = 0, the claim follows immediately from the unitarity of the time 
evolution. In order to prove the induction step, we note that in (A.2), the order of 
differentiation of the wave function on the right hand side is smaller than that on the 
left hand side at least by one. In the case a = 0 and 6 > 0, the induction hypothesis 
yields the inequality 
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and integrating this inequality from 0 to t gives the result. In the case a > 0 and 6 > 0, 
we apply (A.l) together with the induction hypothesis to obtain 

6-1 




< 




bC{l + \tt^) 11 ^ 1 ^- 16=01 




+ cC 


p=0 

6 b 


1 + 


Again integrating over t gives the result. 


□ 
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